Effect of magnetic field on electron neutrino sphere in 

pulsars 



Ashok Goyal Q 

Department of Physics and Astrophysics 
University of Delhi , Delhi - 110 007 India 
Inter University Centre for Astronomy and Astrophysics 
Ganeshkhind, Pune - 411 007 India 

Abstract 

We study the neutrino interaction rates through charged as well as 
neutral current weak interactions in hot, dense magnetized matter. At 
densities near the neutrino sphere, matter in the presence of intense 
magnetic field is polarized and electrons and protons typically occupy 
the lowest Landau level. The weak interaction rates in such a situation 
are severally modified leading to a sizable change in neutrino opacity 
and consequently on the location of the neutrino sphere which also 
gets distorted. This has implications for pulsar velocity explanations. 
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Recently there has been a lot of interest |IJ] - H in the study of neutrino 
transport in dense, hot nuclear matter in the presence of strong magnetic 
fields reported to be present in young pulsars in connection with a possible 
mechanism to explain the recoil velocities of pulsars. Specifically, an asym- 
metry arising due to parity violation in a strong magnetic field can lead to 
asymmetry in the explosion and recoil of newly formed neutron stars. The 
asymmetry could for example, arise in the standard weak interactions of 
neutrinos [|l|] - @, due to neutrino magnetic moment || or due to matter 
induced neutrino oscillations in the presence of magnetic field @] - f|. In 
order to explain the large velocities of neutron stars one needs an asymmetry 
of roughly one percent in the radiated neutrinos which could arise due to a 
number of neutrino scattering reactions on polarised electrons and hadron 
of the magnetized medium. A particularly attractive mechanism to explain 
pulsar kicks is suggested by Kusenko and Segre Q and makes use of resonant 
MSW conversion. The way this mechanism works is as follows: 

The high densities and temperatures reached in proto-neutron star stage 
results in the trapping of neutrinos within their neutrino-spheres i.e. the 
surface from which the optical depth of neutrinos become one. The lo- 
cation of this neutrino-sphere is flavor dependent because of the fact that 
electron neutrinos interact both through charged as well as neutral currents, 
whereas the tau and mu neutrinos have only the neutral current interactions. 
This results in the tau neutrino escaping from deeper layers of star (due to 
their large mean free path) where the temperature is high. The electron 
neutrinos trapped between the tau and electron neutrino sphere can then 
oscillate through resonant conversion to tau neutrinos and become able to 
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escape-thereby making the tau neutrino sphere to effectively coincide with 
the resonance surface. The presence of magnetic field modifies the poten- 
tials relevant for describing neutrino conversion and leads to the distortion 
of resonance surface for MSW conversion through polarisation of the mag- 
netised medium. This leads to tau neutrino emission in different directions 
from regions of different temperatures, therefore with different energies and 
thus leading to pulsar kicks. For this mechanism to work magnetic fields of 
strength B > 10 14 G are required and such high magnetic fields in the context 
of supernova collapse have indeed been proposed in the literature @. 

Such strong magnetic fields would affect the neutrino opacities by modi- 
fying the electron neutrino absorption and scattering cross-sections on nucle- 
ons and leptons which determine the electron neutrino mean free path. The 
location of the v e sphere would thus get shifted and distorted due to asym- 
metry in the cross-section. This effect on the dominant absorption reaction 
v e + n — > p + e~ has been calculated in the literature |8j by modifying the 
phase space distribution of the final state electrons only and by leaving the 
matrix elements unchanged. In the presence of magnetic fields of the order 
of 10 16 G, the motion of the degenerate electrons at densities likely to be 
present near the neutrino sphere (Y e p ~ 10 10 — 10 u g/cc)is in-fact quantised 
and the electrons typically occupy the lowest Landau energy level leading 
to total polarisation of electrons opposite to the direction of the magnetic 
field. Because of charge neutrality viz n e = n p , the protons too are forced 
to occupy the lowest Landau level. In this situation , the matrix elements 
for the absorption and scattering processes from nucleons and leptons get 
modified and have to be calculated by using the exact wave functions for 
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electrons and protons by solving the Dirac equation in the magnetic field 
namely for the lowest Landau level v — 0. Further in order to make numer- 
ical estimates of the neutrino mean free path in hot magnetised matter, we 
require to know the composition of nuclear matter. This is done by consid- 
ering electrically neutral hot n-p-e matter in j3- equilibrium in the presence 
of magnetic field. The nucleons here are treated non-relativistically and at 
the densities of interest, strong interaction effects can be safely ignored. 
The processes that contribute to u e opacity are the absorption 



and the scattering 



+ n(p 2 ) -> p(p 3 ) + e (p 4 ) (1) 

v(pi)+n(p 2 ) -> n(p 3 )+u(p 4 ) (2) 

v(pi)+p(p 2 ) -> p{Pz) + v(Pi) (3) 
v(pi) + e-(p 2 ) -> e~(p 3 ) +v{p 4 ) (4) 

processes. 

The cross-section per unit volume of matter or the inverse mean free path 
is given by 

^ = A- 1 ^) = ^-U n=2 , 3A dp n W fi f 2 (E 2 )(l - f 3 (E 3 ))(l - UE,)) (5) 

where dp n = ^f^E ls ^ ne density of states of particles with four momenta 
p n , f n (E n ) are the particle distribution functions which in thermal equilib- 
rium are given by the usual Fermi-Dirac distributions f n (E n ) = 1/(1 + 
exp{ En ~^ n )) and (1 — f n (E n )) accounts for the Pauli Blocking factor for 
the final state particles. The transition rate Wfi is 

W fl = (2nr5\P f -P l )\M\ 2 (6) 



where \M\ 2 is the squared matrix element summed over initial and final spins. 
In the presence of magnetic field, the density of states are replaced by 

dp nz dp. 



v,0) 



(2vr) 2 2£ n 
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where the sum over the Landau levels is to be performed. If the magnetic 
field is not strong enough to force particles in the lowest v — Landau state, 
the matrix element remains essentially unchanged || and one has to only 
take the proper phase space into consideration. In the presence of quantising 
magnetic field, a situation likely to be present for strong fields at densities and 
temperatures under consideration, the electrons occupy the lowest Landau 
level and the matrix elements have to be calculated by using the exact wave 
functions of relativistic electrons and protons in a magnetic field. In a gauge 
in which the vector potential A = (0, xB, 0), the quantum states are specified 
by the quantum numbers p y ,p z , vands and the energy is given by JO 



E 



m? + p 2 z + 2veB 



(8) 



and the positive energy electron spinor in v = state is given by 

f \ 

E e + m e 




(9) 



electrons being alligned antiparallel to the direction of the magnetic field 
which is along the Z axis. The x-component of the momentum is not con- 
served now and the transition rate becomes 



W 



fi 



(2n) 3 5(E f - E t )6(P fz - P lz )5{Pfy ~ Pry) 



(10) 



V = L x L y Lz is the normalization volume. Squared matrix element summed 
over initial and final spins for processes (l)-(4) in the Standard model is given 
by 

\M\ 2 = 32G 2 F (C v +CA) 2 Pi.P2P3.p±+(Cv-CA) 2 P2.P4Pi.P3-m2m 3 (C^-C 2 A ) Pl .p 4 

(11) 

the vector and the axial vector couplings are given by 



Process Cv C/ 



vn —> pe g v = 1 g& pa 1.22 

vn -> nv -gv/2 ~W 2 

up —y pv 1/2 — 2sin 2 6 w m Pa/2 

ve^ev 1/2 + 2sin 2 6 w 1/2 



During the Kelvin-Helmholtz cooling phase of the proto-neutron star, 
the density is much less than the nuclear density. The nucleons are non- 
relativistic and non-degenerate, the electrons remain relativistic. The neu- 
trinos are non-degenerate and the elastic approximation in the scattering 
processes is quite reliable. The phase space integrals in eqn(5) can then be 
carried out analytically ]TT and we obtain the mean free paths as 

A -< 0) = V^ + ^^ i + ^Vq-/^ (12) 

K l (0) = ^{g 2 v + ig\)El„ N (13) 

= -^9lE 2 n p (14) 
or 12 

K\0) = ^f(C 2 v + C 2 MTE 2 u (15) 



Considering now the magnetic field that affects only the electrons and 
which is not strong enough to confine the electrons to the lowest Landau 
state, the matrix elements for the processes remain essentially unchanged 
and modifying only the phase space, the important absorption cross-section 
turns out to be |§ 



Cabs(B) G 2 F , 2 2 . 1 
^ 9v + 39A)eB l + exp(E u + Q-^)P 

(E v + Q) y £(2-6 u ,o) ' 



V 



{E v + Q) 2 - m? e - 2veB 



(16) 



The cross-section for electron neutrino scattering in hot dense matter in 
the presence of magnetic field has been derived in 0. The other processes 
remain unaffected by the magnetic field. 

In the interesting case of strongly quantised magnetic field when both 
electrons and protons in stellar matter are polarised with their spins aligned 
anti parallel and parallel to the direction of the magnetic field respectively, 
the matrix elements for the above processes are calculated using exact wave 
functions for v — Landau state .For non-relativistic nucleons and relativistic 
leptons, we get for the matrix element squared 



M 2 ahs 



8G 2 F m 2 m 3 (E 4: + p Az 



(9v + 9 a? + + (g v - g A f + \g\E x CosQ 



exp 



2eB 



(Plx + Plxf + (P3y +P4yf 



M. 



■pv 



16G 2 F m 2 2 (pi.p4 + 2p 1Z -P4z) 



1 



Mi 



exp 



IQG 2 F 



2eB 

(C 2 V + C\ 



{Pix - Plxf + {Piy - Ply) 2 

(E X E A + pizP2z)(E 2 E 3 + p 2 zP4z) 



(17) 



(18) 
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-(Emz + E 4 p lz )(E 2 p 3Z + E 3 p 2Z ) 



(E ± E 4 + PizP4z)(E 2 P3z + E 3 p 2Z ) 



-{Eip±z + E 4 p lz )(E 2 E 3 + p 2Z p 3Z ) 
1 



exp 



2eB 



(PAX-Plx) +(P4y-piyY 



The scattering cross- sect ions are now given by 



Cscatt(B) 

V 



■^scatt 
1 



d 3 p 4 



2E 1 L X J (2tt) 3 2£W (2n) 2 2E 2 
(2nf5(P Y )5(P z )5(E)\M\ 2 

f 2 (E 2 )(l - f 3 (E 3 ))(l - f 4 (E 4 )) 



(19) 



dp 2Y dp 2Z f dp 3Y dp 3Z 



(2ir) 2 2E 3 



(20) 



and the absorption cross-section is obtained by exchanging the fourth and 
the second particles. The integrals over y-component of the momenta can be 
done by using the y-component momentum conserving delta function. The 
rest of the integrals can be done by using 5(E) and S(P Z ) in the limit of 
non-degenerate, non-relativistic nucleons. The exponential functions in the 
squared matrix elements can be approximated by one in this limit and further 
in the limit of elastic scattering approximation. Carrying out the integrals 
we obtain 



— y — = Kbs{B) 



An 



eBn 



1 



1 + exp{-(E v + Q- fi e )/3} 



{{9v + 9a? + 4g 2 A } + {(g v + g A f - ±g\}cos6 

(21) 
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cr vp (B) 
V 



2G 2 F C\eB r°° 



(27T) 



(J a en r°° 

d P2Z f 2 (E 2 )(E u 



Viz 
2m p 

|2 _ \ I J? _|_ P2Z \-j _ PnwZ + 

3m p " 2m p m p 



gig 
2m„ 



which can be approximated by 



V 



TV 



n p E u 



where 



2eB r°° 



n 



p (2tt) 2 



/ f2(E 2 )dp 2 z 
Jo 



Vve(B) AGpeB 
V ~ (2tt) 3 



(C£ + C^)(l + cos 2 9) - AC v C A cos9 



1 + exp(-pL e /T) 



(22) 



(23) 



(24) 



(25) 



The neutrino-neutron scattering cross-section remains unchanged except 
through its dependence on neutron density which gets affected by the mag- 
netic field. 

For the noninteracting, hot n-p-e matter in beta-equilibrium we have the 
number densities given by 



2eB 



(27T) 2 
-(/i e -> -fi e ) 

2eB 



£(2 - S v , ) / dp z - 

} . JO 1 



f 



+ exp{f3(E ev - /j, e )} 



(26) 



n r . 



n n = 



(2vr) 2 

— / p^ dp 

n 2 Jo 1 + exp{/3(E nu - /i n )} 



+ exp{[3(E pu - fi p )} 



(27) 
(28) 



The energy expressions in the presence of magnetic field are given in 
eqn(9). Charge neutrality requires n e = n p and /^-equilibrium relates the 
chemical potentials viz. fi n = fi p + /i e 

In Fig.l we have plotted the proton fraction Y p = - n ^_ n as a function 
of density at two different temperatures 1 MeV and 10 MeV respectively 
for different values of the magnetic field. The neutron fraction is given by 
Y n = 1 — Y p . At T=l MeV, there are very small number of positrons present 
and the effect of the magnetic field is to raise the proton fraction so much so 
that at B = 10 4 MeV, the matter becomes predominantly proton rich. This 
had the effect of not only modifying neutrino cross-sections quantitatively 
but to change the relative importance of various processes in the presence of 
magnetic field. In Fig. 2 we have plotted the neutrino mean free path in meters 
as a function of baryon density for different values of the magnetic fields for 
matter at T=l MeV. The neutrino energy is taken to be E\ = 3T. We observe 
as expected that neutrino absorption cross section dominates the mean free 
path even in the presence of highest magnetic fields considered here as it 
does in the free field case. We also find that the absorption mean free path 
rises rapidly with density since at low temperatures due to Pauli-Blocking 
effect, the absorption process gets suppressed. At higher temperatures, there 
are lots of positrons present and electron degeneracy decreases and the effect 
is absent as can be seen from Figures 3 and 4 where we have plotted the 
mean free paths for matter at T=10 and T=30 MeV respectively. In fact 
at these temperatures, the free path decreases with increase in density. The 
absorption and neutrino electron scattering cross sections develop asymmetry 
when the matter is polarised as can be seen from eqns.(21) and(25). The 
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asymmetry in the absorption cross section for neutrinos propagating along 
and opposite to the direction of the magnetic field is roughly of the order of 
20 %. The asymmetry is far more pronounced in neutrino electron scattering 
, so much so that scattering cross section for neutrinos in the direction of the 
magnetic field is ten times more compared to neutrinos progating parallel to 
the field. Thus at low temperatures, when the matter is highly polarised the 
neutrino sphere itself gets severely modified and distorted. This results in 
neutrinos along the magnetic field being emitted from deeper layers where 
the temperature is higher in comparison to neutrinos which are emitted anti 
parallel to the field direction. A change in cross sections ( by an order of 
magnitude even at high temperatures ) caused by a large magnetic field 
will change the density near the neutrino sphere considerably. We can see 
from Figures 2, 3 and 4 that increase in neutrino cross section is roughly 
proportional to B ( albeit for polarised matter ) thereby shifting the neutrino 
sphere to larger radii and hence to lower densities. This will have implications 
for the implementation of Kusenko-Segre mechanism to explain the observed 
pulsar velocities. 
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Figure Captions 

Fig.l Proton fraction Y p as a function of baryon density in gm/cc for 
matter at T = 1 MeV and T = 10 MeV respectively The solid, long dashed, 
dashed and dotted curves are for B = 0, 10 2 , 10 3 and 10 4 MeV 2 respectively. 

Fig. 2 Neutrino mean free path in meters as a function of baryon density 
for matter at T=l MeV and B = 0, 10 2 and !0 4 MeV 2 . The neutrino energy 
is taken to be 3T. The solid curve is for the absorption process and the long 
dashed, dashed and dotted curves are for neutrino scattering on neutrons, 
protons and electrons respectively. 

Fig. 3 Neutrino mean free path in meters as a function of baryon density 
for matter at T=10 MeV and B = and 10 4 MeV 2 . The neutrino energy is 
taken to be 3T. The solid, long dashed, dashed and dotted curves are as in 
Figure2. 

Fig. 4 Neutrino mean free path in meters as a function of baryon density 
for matter at T=30 MeV and B = and 10 4 MeV 2 . The neutrino energy is 
taken to be 3T. The solid curve, long dashed, dashed and dotted curves are 
as in Figure2. 
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Figure 1: 
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Figure 3: 
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Figure 4: 
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